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V  IntroducliQn 

'"‘Langmuir-Blodgett  film  monolayers  are  formed  when 
amphiphilic  molecules  in  solution  are  deposited  on  water 
surface  followed  by  compression.  During  the 
compression,  the  randomly  oriented  molecules  reorient 
themselves  with  the  hydrophilic  head  groups  on  the 
water  surface  and  the  hydrophobic  tail  standing  away 
from  the  water  surface.  The  organic  ultra-thin  films 
with  highly  ordered  layer  structure,  have  shown  a  wide 
range  of  potential  applications  in  biochemistry,' 
microelectronics,  integrated  optics,  microlithography, 
and  many  other  fields  It  is  very  Important  to  get  a 
better  understanding  of  the  relationship  between  the 
properties  of  the  films  and  the  micro-structure,  for 
example,  the  molecular  conformation,  packing, 
orientation  and  motion  in  the  films.  For  this  purpose, 
within  the  past  few  years,  different  computer  simulation 
methods  yvere  applied  on  LJjP'film  monolayers,  including 
Monte  Carlo  niethods  ^5]  and  molecular  dynamics 
simulations  In  this  paper,  we  present  the  results 

of  a  molectiUfr  dynamics  simulation  of  a  fatty  acid 
■  (CH3(CH2)27COOH)  monolayer  on  a  water  surface.  This 
fatty  acid  is  a  preliminary  model  for  our  next  study 
object:  CH3(CH2) ,  5C«C-C«C(CH2)8COOH. 

Methodology 

QUANTA  2.0  (CHARMm  21),  a  comprehensive, 
Integrated  molecular  modelling  package  developed  by 
Polygen  Corporation  is  used  in  the  present  simulation.  In 
this  calculation,  the  Steepest  Descents  method  (7)  is 
employed  for  energy  minimization  and  the  Verlet 
algorithm  is  employed  for  the  dynamics  simulation.  The 
hardware  employed  is  a  GS1000  Graphics  Supercomputer 
from  Stellar  Computer  Inc. 

The  simulation  starts  with  169  molecules,  in  a 
13x13  matrix,  parallel  to  each  other  standing  vertically 
on  a  layer  of  water.  During  the  dynamics  simulation,  the 
position  of  the  water  molecules  are  frozen,  A  15x15 
molecule  matrix  (225  molecules)  without  water  is  also 
studied. 

To  relax  the  conformation  and  remove  bad  contacts, 
the  system  undergoes  energy  minimization  at  zero  K.  The 
potential  energy  terms  concerned  are: 

"  ^bond  length  +  ^bond  angle  +  ^dihedral  angle 

+  ^Improper  torsion  +  ^electrostatic  +  ^van  der  waata 

The  energy  for  hydrogen  bonding  is  included  in 
electrostatic  and  van  der  waal  terms. 

The  optimized  system  is  then  heated  to  the  desired 
temperature.  I.e.  300K,  by  gradually  assigning  random 
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velocities  to  each  atom.  The  relationship  between  {he 
atom  velocity,  the  kinetic  energy  and  temperature  of  the 
system  is: 

Bkinoilo  -  £  m,(v2}  .  3.  NkoT 

2  1-1  2 

where  N  is  the  total  number  of  atoms  and  m,  and  v,  are 
the  ith  atom's  mass  and  velocity,  respectively;  kB  is 
Boltzmann’s  constant. 

When  the  system  is  heated  to  300K,  It  undergoes  a 
equilibration  period  to  stabilize  average  temperature  and 
structure.  After  the  system  is  equilibrated,  the  final 
molecular  dynamics  simulation  starts.  The  coordinates 
of  the  molecules  are  recorded  periodically.  The  sets  of 
coordinates  recorded  during  the  dynamics  simulation  are 
then  analyzed  to  extract  information  about  molecular 
packing,  orientation  and  motion. 

Results  and  Discussion 

Water  strongly  limits  the  mobility  of  the  molecules, 
especially  the  polar  head  groups.  In  Table  I,  the  vertical 
mobility  of  the  two  end-groups  of  the  molecules  Is 
represented  by  the  fluctuations  in  height.  Without  water, 
both  end-groups  have  similar  fluctuations  as  shown  in 
Table  I  and  Fig.  la.  When  water  is  introduced  into  the 
system,  the  surface  formed  by  the  polar  groups  (the 
lower  curve  in  Fig.  1b)  is  more  even  due  to  their 
hydrophilic  interaction. 

During  the  simulation  the  molecules  tend  to  adjust 
the  Intermolecular  distance  to  reach  the  minimum 
potential  energy.  For  closely  packed  systems  (i.e.  with 
smaller  area  per  molecule),  the  final  molecular  areas  of 
21 -22 A2  (see  Table  I)  are  consistent  with  that  for 
'close-packed  heads’  phase  of  fatty  acids  monolayers  [8]. 
For  closely  packed  system,  there  Is  pseudohexagonal 
packing  with  defects  (Fig.  2).  Relatively  loosely  packed 
systems  however  have  significantly  less  packing  order. 

Comparing  Fig.  3a  and  3b  we  can  see  that  the  center 
molecule,  surrounded  by  other  molecules,  has  low 
mobility  and  flexibility,  represented  by  the  tilt  and 
molecular  length  respectively;  the  molecule  at  a  corner 
of  the  molecule  matrix  has  more  mobility  and  flexibility. 
When  loosely  packed,  there  is  more  space  around  each 
molecule  and  the  cohesion  is  weaker,  therefore  the 
molecules  will  have  more  tilt  (see  Table  I)  and  stronger 
edge  effect  as  shown  in  Fig.  4:  the  tilts  for  central 
molecules  are  less  then  10"  and  for  molecules  at  the 
matrix  corner  can  be  as  high  as  60”.  To  correct  for  the 
edge  effect,  simulations  with  periodic  boundaries  are 
being  made  and  will  be  discussed. 

1  Sanslualan 

Simulations  of  the  LB  films  without  periodic 
boundaries  show  clearly  the  effect  of  water  and  packing 
condition  on  the  behavior  of  molecules  in  the  monolayer. 
If  is  instructive  to  see  appearance  of  domain  boundaries 
and  lattice  defects  in  the  simulation  process  as  result 
already  seen  experimentally  but  whose  origin  Is  not  well 
,  understand.  The  strong  edge  effect  makes  it  necessary  to 
lintroduce  periodic  boundaries  in  future  simulations. 
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Fig.  2  The  projection  of  average  mass  centers  of  the 
molecules  (open  circle)  on  a  water  surface  with 
molecular  area  22A2/molecule  (starting  from  25A2). 


Table  I.  The  height  fluctuation  [Aj  of  the  end-groups,  the 
change  of  molecular  area  (A2/molecule]  and  the  average 
tilt  [°]  in  the  dynamics  simulations  with  and  without 
water.  The  simulation  time  scale  is  pico  second  [ps] 

Simulation  15x15  molecules  13x13  molecules 
without  water  with  water 

Molecular  Area  [A2/molecule] 


Start 
End 

Head  Group 
Tail  Group 

First  30ps 
Last  30ps 


16  36  25 

21  25  22 

Fluctuation  In  Height  [A] 

2.52  4.18  1.14 

2.93  4.16  3.83 

Average  Tilt  [°] 

9.37  14.65(1 7ps) 


36 

27 

0.86 

3.54 


12.35  20.02 

18.44  30.00 
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Fig.  3  The  molecular  length  [Aj  and  tilt  [”]  as  a  function 
of  dynamics  sinvlation  time(in  pico  seconds)  for  (a)  a 
center  molecule  and  (b)  a  corner  molecule  In  a  loosely 
packed  system(27A2/molecule).(  0*  tilt  means  the 
molecule  is  vertical  to  water  surface). 
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Fig.  1  The  time  averaged  height  of  the  mass  center  and 
the  two  end-groups  of  the  molecules,  (a)  Without  water 
(21  A*/niolecule).  (b)  With  water  (22A2/molecule). 


Fig  4  The  average  tilt  angle  (»]  of  molecules  in  a  loosely 
packed  system  vs.  the  molecule  position  in  two 
dimensions  with  molecular  area  27A2  for  30  os 
simulation. 
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